Detecting atmospheric neutrino oscillations in the ATLAS detector at CERN 
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We discuss the possibility to study oscillations of atmospheric neutrinos in the ATLAS experiment 
at CERN. Due to the large total detector mass, a significant number of events is expected, and 
during the shutdown phases of the LHC, reconstruction of these events will be possible with very 
good energy and angular resolutions, and with charge identification. We argue that 500 live days of 
neutrino running could be achieved, and that a total of ~ 160 contained Vy, events and ~ 360 upward 
going muons could be collected during this time. Despite the low statistics, the excellent detector 
resolution will allow for an unambiguous confirmation of atmospheric neutrino oscillations and for 
measurements of the leading oscillation parameters. Though our detailed simulations show that 
the sensitivity of ATLAS is worse than that of dedicated neutrino experiments, we demonstrate that 
more sophisticated detectors, e.g. at the ILC, could be highly competitive with upcoming superbeam 
experiments, and might even give indications for the mass hierarchy and for the value of 6n. 

PACS numbers: 14.60.Pq,13.15.+g 



INTRODUCTION 



Neutrino oscillation physics is entering the era of pre- 
cision measurements, and a plethora of new dedicated 
experiments are being designed, or already under con- 
struction. The most widely discussed technologies are 
reactor experiments 0, , superbeams [1, Q| , advanced 
atmospheric neutrino detectors 0, @] , beta beams 0, [1] , 
and neutrino factories [T3| ■ In this article, we enter- 
tain the possibility to use the ATLAS detector at CERN 
to study oscillations of atmospheric neutrinos, an idea 
which has first been brought up by F. Vannucci [Til. H3|. 
A significant number of atmospheric neutrino interac- 
tions will take place in the 4 kt hadronic calorimeter of 
ATLAS, and can be reconstructed with excellent energy 
and angular resolution during phases where the LHC is 
not running at full luminosity. Moreover, the magnetic 
field will allow for the discrimination between neutrinos 
and anti-neutrinos. 

To study the physics that can be done with these event, 
we have developped a code that allows for the simulation 
of event spectra, taking into account a full three-flavor 
treatment of neutrino oscillations, the finite energy and 
angular resolutions of the experiment, and the detector 
geometry. To derive high-level results, such as parameter 
sensitivities, from the simulated data, we use a x 2 anal- 
ysis including systematical uncertainties, parameter cor- 
relations, and degeneracies. 

The paper is organized as follows: After briefly review- 
ing the physics of atmospheric neutrinos in Sec. [2l we 
discuss the capabilities of ATLAS to detect atmospheric 
neutrinos in Sec. [3] We then describe the technical de- 



tails of our simulations in Sec. [H and present the results 
for the sensitivity to the leading atmospheric oscillation 
parameters, to the mass hierarchy, and to other three- 
flavor effects, in Sec. [5] Our conclusions will be presented 
in Sec. [U 



2. OSCILLATIONS OF ATMOSPHERIC 
NEUTRINOS 

Atmospheric neutrinos are produced by interactions of 
cosmic rays with the atmosphere at a height of around 10 
to 20 km above ground. The main production reaction 
is 



p + Atomic Nucleus 



Further Hadrons 



(1) 



e + + v e 



Due to strong 



oscillations, the flux of up- 



ward going neutrinos is significantly lower than that 
of downward going neutrinos. This leading oscilla- 
tion effect has been unambiguously detected experimen- 
tally [n, lh, m, m, ■ 

In contrast, subleading three- 



flavor effects, such as u u 



oscillations driven by 6*13 
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and Amj! , are very hard to detect with atmospheric neu- 
trinos because the flux ratio (y^ + P^)/(t / e + v e ) ~ 2 (cf. 
Eq. |l|)), in combination with the close-to-maximal value 
of 623 leads to a cancellation of — > v e and v e — > 
oscillations (l8j . 

Atmospheric neutrinos cover a wide range of energies, 
but their flux decreases rapidly for E > 3 GeV. The dis- 
tances travelled by them before detection range up to 
12,742 km, the diameter of the Earth, for upward going 
neutrinos. Detailed calculations of the atmospheric neu- 
trino fluxes have been performed by Honda et al. [l9l ]. 
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Battistoni et al. [H, and Barr et al. [22J. In our 
study, we will use the Honda fluxes [l9[ . 

ATLAS will be sensitive mainly to the leading — + v T 
oscillations, but due to the difficulties associated with the 
reconstruction of v T interactions, the most important os- 
cillation channel will be — > disappearance. Neglect- 
ing sublcading effects driven by Am|i, the corresponding 
oscillation probability in matter is given by [23| : 

P w = 1 - sin 4 9 23 sin 2 26-™ sin 2 Ci 3 A 

- - sin 2 26»23 [l - cos(l + A) A cos C 13 A 

+ cos26»^sin(l + A)AsinC*i 3 A] (2) 

Here, we have used the notation from (23[: A = 
Am^L/AE is the oscillation phase, A = 2EV/Am 31 de- 
scribes matter effects, sin 26^ = C 1 - 3 1 sin26'i3 is the mix- 
ing angle in matter, and C13 is given by C13 = [sin 2 28i 3 + 
( cos 26»i3 ~ A) 2 ] 1 / 2 . The first line in Eq. © describes the 
disappearance due to matter-enhanced — > v e oscil- 
lations, while the other terms correspond to — > v T 
oscillations. A generic three-flavor effect is the strong in- 
fluence of matter on the v u — > v., and v u — * v T channels 
in some regions of the parameter space [2J, |25j . In the 
vacuum case, i.e. for A = 0, C13 = 1, and 6^3 = 6*13, 
Eq, ([2]) reduces to 

P WI = 1- [sin 4 6» 23 sin 2 26>i 3 -cos 2 13 sin 2 26> 23 ] sin 2 A, 

(3) 

i.e. to two-flavor oscillations with a mixing angle close to 
#23. For the neutrino oscillation parameters, we adopt 
the following best-fit values in our simulations (see e.g. 
[M 11, 11): 



sin 2 26»i2 


= 0.79, 


sin 2 20 13 


= 0.12, 


sin 2 26*23 


= 1.0, 


<5cp 


= 0.0, 


Am 2 -! 


= 8.1 • 10~ 5 cV 2 


ZW31 


= 2.2 • 10~ 3 eV 2 



(4) 



3. RECONSTRUCTION OF NEUTRINO 
INTERACTIONS IN ATLAS 

The ATLAS detector at CERN has an onion shell 
structure, similar to that of most other modern col- 
lider experiments: The inner high-resolution tracking 
detectors are surrounded by the electromagnetic and 
hadronic calorimeters, large superconducting magnets, 
and a muon tracking system [3(|. ATLAS has a total 
mass of 7 kt, but part of it is attributed to the non- 
active support structure of the experiment. Reconstruc- 
tion of neutrino interactions will only be possible if the 
energy and direction of the secondary charged lepton can 



be seen, and for a good energy and angular resolution, it 
is desirable to also reconstruct the energies of the other 
interaction products. For interactions, this requires 
the interaction to take place inside the 4 kt hadronic 
calorimeter, so that the muons will travel through the 
whole muon system. Furthermore, the neutrino must 
have sufficient energy to yield a sizeable signal. We esti- 
mate that a threshold of 1.5 GeV should be realistic. For 
v e and u T , the track of the secondary charged lepton can 
only be seen with good resolution if the interaction takes 
place in the inner tracking detectors, which, however, 
have negligible mass. Energy reconstruction for v e and 
v r should in principle be possible also for interactions 
taking place further outside, but the resolution would 
be quite poor since only information from the hadronic 
calorimeter could be used. Therefore, we do not take 
v e and v r into account at all. We do, however, include 
so-called upward going muons, i.e. muons created by neu- 
trino interactions in the rock below the detector. More 
precisely, we consider muons coming from the zenith an- 
gle range —1 < cos6> < —0.1 to be induced by neutrinos. 
Note that ATLAS, being a magnetized detector, is able 
to determine the muon charge and can thus discriminate 
between neutrinos and anti-neutrinos. 

A very important issue for the detection of atmo- 
spheric neutrinos in ATLAS is triggering. The charac- 
teristic signature of an atmospheric neutrino is a high- 
energy charged lepton track originating from a vertex 
inside the detector, without any visible ingoing particle. 
During normal LHC operation at full luminosity, such 
signals will be cloaked by the pile-up of pp interaction 
products, and even if they could be seen, they might 
still be confused with decay products of neutral hadrons. 
Besides, the ATLAS trigger will only be sensitive during 
the bunch crossings, because the time between them is 
required to read out the accumulated data. For these 
reasons, ATLAS can be used for neutrino physics only 
during the ~ 200 days per year where the LHC is not 
running in collider mode [3 if . This number arises from 
the 14 week winter shutdown, and from regular short 
maintenance shutdowns throughout the year. Moreover, 
even during LHC operation, the pp interaction rate in 
the detector will be sufficiently low during the ramp- up 
and ramp-down phases. Of course, also the ATLAS de- 
tector will require maintenance, so that only part of the 
aforementioned time windows will be available for neu- 
trino physics. If ~ 100 days of neutrino running per year 
are feasible, this would yield ~ 160 contained events and 
~ 360 upward going muons within 5 years. We will show 
in Sec. [5] that this is sufficient to detect neutrino oscil- 
lations and to perform measurements of the oscillation 
parameters. 

Note that the number of neutrino interactions expected 
in CMS should be even larger than that in ATLAS, since 
CMS has a total mass of 12 kt. [3^|. However, most of this 
mass is concentrated in the massive iron return yokes of 
the magnet, so that all non-muonic interaction products 
will be scattered several times before reaching the active 
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detector components. This will inhibit a reliable recon- 
struction of the primary neutrino properties, and hence 
we will not consider CMS further in this article. 

Since the leading oscillation channel for atmospheric 
neutrinos is — > v T , and the r has a 17% chance of 
decaying into fxv^v T |33j, one has to ask the question 
whether v T interactions, followed by a leptonic r decay, 
can be clearly separated from simple charged current 
interactions. However, this should be possible, since typ- 
ically, the muon from the r decay will have a relatively 
low energy, which is most probably below the threshold 
used in the analysis. 



4. SIMULATION OF ATMOSPHERIC 
NEUTRINOS 

To study the prospects of doing neutrino oscillation 
physics with ATLAS, wc have developed a simulation 
code for the calculation of event rates as a function 
of the neutrino energy and zenith angle, and for the 
X 2 analysis of the simulated data. We obtain the 
binned event rates by folding the initial neutrinos fluxes 
$// (E*, L(6 l v )) , [H], the three-flavor oscillation proba- 
bilities in matter P(f — > /, E*, L{6 l v ), 9), the cross sec- 
tions <jf(E*) H3, [35[, and a detector response function 
R^ (E„, 6 V ) according to the formula 

N} j =Mj2R ij (E^el)-a f (E^ 

k,i 

In this expression, E* and 6 V denote the neutrino en- 
ergy and zenith angle at the sampling point with in- 
dices (fc,/), while stands for the binning used in 
the analysis. The number of bins is roughly chosen 
according to the energy and angular resolutions (see 
Tab. [H. 1 L{0 l v ) is the distance travelled by neutrinos 
coming from the direction 0\,, and represents the vec- 
tor of oscillation paramaters (including the matter poten- 

tial), i.e. = (0i2,0 1 3,023,<fcp,Am2 1) Am§ lJ V). The 
cross sections used in our code cover an energy range 
from 100 MeV to 1 TeV. For neutrinos with higher ener- 
gies (up to 10 TeV), which can contribute only to the 
upward going muon sample, we extrapolate the cross 
sections by making the assumption that a(E)/E is con- 
stant at such high energies. We only use total charged 
current cross sections. The detector response function 
R^(E^,9l) determines which fraction of neutrinos with 



1 This binning could not be employed in the analysis of real data, 
since the event numbers in each bin would be very small, but for 
our purposes, it gives a good estimate of the information that is 
contained in the data. 



initial energy E„ and zenith angle 9 l v is reconstructed into 
bin Since the detector response to neutrinos has 

not been studied in full detector Monte Carlo simulations 
yet, we resort to approximating R by a double Gaussian 
resolution function, multiplied with an efficiency factor 

e(E u ,e u y. 

■z-Jo ^ exp l" M(^) J" 

Here, <je{E v ) and a a (E u ) denote the energy and angu- 
lar resolutions, respectively. a{6 3 r , 4> r , l v ) is the angle be- 
tween the initial neutrino direction, parameterized by the 
zenith angle 8^, and the reconstructed direction, which 
is defined by the zenith and azimuthal angles 9 3 r and (f> r ; 
Ze and Z a are normalization factors which ensure that 
the total number of events is conserved in the "smearing" 
of the spectrum. Due to lack of Monte Carlo results on 
the reconstruction efficiency in ATLAS, we omit e(E v , 9 V ) 
for contained events. For upward going muons, we use e 
to scale the flux according to the target volume, which 
varies with energy and zenith angle (see appendix [A] for 
details). 

For the resolutions, we use estimates based on the ex- 
perience of the Super-Kamiokande collaboration [l5l . l36ll 
and on the information from the ATLAS proposal [301 ]. 
In Super-Kamiokande, where only the secondary muon 
is seen, the energy resolution for contained Multi-GeV 
events is 17%, and the angular resolution is 17° above 
1.5 GeV. The ATLAS calorimeter measures the full de- 
posited energy with a resolution of order <je/E ~ 
50%/y/E/GeV, while the muon system has an energy 
resolution on the per cent level. The angular resolution of 
the muon system is excellent because track resconstruc- 
tion is possible, while the directional information from 
the hadronic calorimeter is severely limited by its rough 
segmentation. We will assume that the excellent muon 
reconstruction, in conjuction with the information from 
the hadronic calorimeter, will yield an overall neutrino 
energy resolution of 5%, and an angular resolution of 7° 
for contained events. Let us emphasize that these num- 
bers are only estimates, and for more reliable numbers, 
detailed detector Monte Carlo simulations are indispens- 
able. For upward going muons, no energy reconstruction 
is possible, and the angular resolution is essentially given 
by the average difference between the neutrino and muon 
directions. Based on [l5| . we take a a — 5° for upward 
going muon events. 

To account for the case that our estimates are not con- 
servative enough, we will also consider a scenario where 
the resolutions are as poor as in Super-Kamiokande, and, 
in addition, the live time is taken to be only 250 days. 
On the other hand, we will also discuss a very opti- 
mistic setup with oeIE = 2% and <r a = 2° for contained 
events. In this setup, we also lower the energy threshold 
to 0.3 GeV and increase the exposure time to 2,000 days. 
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Ixunning time 
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9^0 rl 
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Energy threshold 


1.5 Gev 


1.5 Gev 


(J.o Gev 


&e / E (contained events) 


r.C/ 


1 ( /o 


9^ 


o~ a (contained events) 


r-O 
i 


1 1° 
1 f 


9° 


a a (upward muons) 


5° 


5° 


5° 


E bins 


30 


20 


90 


cos 6 bins (contained events) 


36 


10 


80 


cos 6 bins (upward muons) 


27 


27 


27 



Table I: Properties of the realistic and conservative ATLAS 
scenarios, andof an assumed ILC detector. 



These assumptions will most probably not apply to AT- 
LAS, since they would correspond to a situation where 
neutrino reconstruction is possible even during normal 
LHC operation, and the reconstruction efficiency is al- 
most 100%. However, they may be interesting in the 
context of future projects such as the proposed ILC de- 
tectors [37j . In an ILC experiment, no pile-up of hadronic 
events will occur, so neutrino detection is not restricted 
to the maintenace phases of the accelerator. The para- 
meters of our three scenarios are summarized in Tab. U 
To analyze the simulated event spectra, we use a \ 2 
fit which distinguishes neutrinos from anti-neutrinos, i.e. 
we assume 100% charge identification efficiency. For each 
species, we define a \ 2 function of the form 

+ 2N ij (Q tIue )hx( Nij ^ true l ) 

vzij(e fit ,3)y 

± Xpull,osc ± Xpull.sys" 

Here, iVtj(0tnie) is the "observed" event rate in bin 
for the "true" oscillation parameters Otrue = 
(#12, 013) 023, S C p, Am? 21 ,Aml 1 , V), and T tJ (<d m , a) is the 
event rate that would be expected for the hypothesized 
parameters Gst, and for biases a arising from systemat- 
ical errors in the experiment. We have introduced pull 
terms 

2 _ (@fc,fit — ©A;, true) 2 /„s 

Xpull,osc / / 2" ' V / 

k °©* 

to account for external input on the oscillation para- 
meters. a® k determines how strongly a fit value far from 
the externally given one is disfavoured. Pull terms are 
provided only for those parameters which are marginal- 
ized over in the fit. We take relative la uncertainties of 
10% for 6>i2 and 23 , 5% for Am|i, and 30% for Am^. 
6>i3 is assigned an absolute uncertainty of 10°. 



Error Type <j a 

Overall normalization for contained events 20% 

Relative normalization for anti-neutrinos 5% 

Normalization for upward going muon events 20% 

Tilt of E spectrum 5% 

Tilt of 6 spectrum for contained events 10% 

Tilt of 9 spectrum for upward going muon events 2% 



Table II: Systematical errors in ATLAS. 

Pull terms for systematical biases are similar to those 
for the oscillation parameters: 

a 2 

2 = V — ^ (71 

Apull.sys 2 ' V ' > 

k a " 

They disfavor fit values afc which are further from zero 
than can be expected from the systematical uncertainties 
<j ak . We have defined a such that the case of vanishing 
systematical errors corresponds to a = 0. The various 
types of systematical errors we consider, are summarized 
in Tab. [TTJ Besides the usual normalization errors we 
allow also for "tilts" in the event spectrum, which are a 
simple way of accounting for energy or angle dependent 
biases. 

Note that many of the technical details of our simu- 
lations and of our statistical analysis procedure follow 
ideas that have been previously realized in the GLoBES 
code (HI, EH- 

5. SIMULATION RESULTS 

In this section, we present our main results on the sen- 
sitivity of ATLAS to the leading atmospheric oscillation 
parameters and to three-flavor effects. In Fig.[TJ we com- 
pare the expected confidence regions in the 023-Am^ 
plane for the three scenarios introduced in Sec. [H and in 
Fig. [2] we relate them to the results of a global three-flavor 
fit to existing atmospheric neutrino data [26|. |40| . and to 
the expected performance of T2K fill, |42u43| . From the 
plots, we first notice that in all three scenarios, neutrino 
oscillations can be confirmed at better than 3<r, i.e. the 
case 023 = 0, Am§i = can be ruled out. However, 
in the conservative case, the sensitivity is only marginal 
due to the limited statistics and resolutions, and could 
be easily spoiled by unfavorable statistical fluctuations 
in the real experiments. 

The precision with which the mixing angle can be de- 
termined is poor in all three scenarios: Even in the op- 
timistic ILC case, the uncertainty at 3<r is about ±30%, 
which is still slightly worse than the current bound from 
Super-Kamiokande, K2K, and MINOS(~ 23%, cf. (Hi). A 
measurement of 023 in ATLAS is not competitive at all. 
The sensitivity to 02 3 in atmospheric neutrinos comes 
mainly from the up-down asymmetry of the neutrino flux, 
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Figure 1: Sensitivity of ATLAS and of a possible ILC detec- 
tor to the leading atmospheric oscillation parameters #23 and 
Am^i. The shaded regions correspond to the la, 2a, and 
3a confidence levels for the realistic scenario, while the dark 
contours correspond to the conservative case. For the ILC 
scenario (light-colored contours), we show only the 3a con- 
tour for clarity. A normal mass hierarchy has been assumed 
in the plot, but we have checked that all results except the 
sensitivity to the mass hierarchy, are similar for the inverted 
hierarchy (see text for details). 
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Figure 2: Sensitivity of ATLAS to the leading atmospheric 
oscillation parameters #23 and Am^ in comparison with a 
global fit to existing atmospheric neutrino data (dominated 
by Super-Kamiokande) [26tl4(il|. and with the expected perfor- 
mance of T2K. For the latter experiment, we show only the 3a 
contour for clarity. A normal mass hierarchy was assumed in 
the plot, but we have checked that the results for the inverted 
hierarchy are analogous. 



which can only be measured accurately with good statis- 
tics. The excellent resolution of ATLAS cannot compen- 
sate its small mass here. 

The situation is better for Am^: For the realistic 
scenario, ATLAS shows a performance similar to Super- 
Kamiokande, at least at low confidence levels. The 3tr 
contour, however, extends much further to large Am 31 . 
In the ILC scenario, the sensitivity to Ato^ is even com- 
parable to that of T2K. An atmospheric neutrino oscil- 
lation experiment is sensitive to the mass squared differ- 
ence mainly through the shape of the zenith angle spec- 
trum at around cos 6 0. These directions correspond 
to baselines below the first oscillation maximum, where 
the depletion of the flux is just setting in. Wash-out due 
to poor detector resolutions can greatly limit this mea- 
surement, so AT LAS-like detectors have a considerable 
advantage here. 

Turning to the discovery potential for generic three- 
flavor effects, we remark, that although a very large value 
of sin 2 2#i3 = 0.12 was assumed in Fig. [1] and Fig. [2j the 
wrong hierarchy solution 0, HEI cannot be fully ruled 
out at 3(7 even by an advanced ILC detector. In the two 
ATLAS scenarios, there is no sensitivity to the hierarchy 
at all. If the true hierarchy is inverted, matter effects 
are shifted to the anti-neutrino channels, where they are 
cloaked due to the smaller anti- neutrino cross sections, 
and we have checked that in this case, even the ILC sce- 
nario has only a la discovery potential. One might have 
hoped that matter effects in the v ^ — > Vn channel could 
have helped to resolve this degeneracy [24], Hfl, but it 
turns out that the event numbers are too small in the 
relevant energy range of several GcV. 

For the same reason, we also expect a poor perfor- 
mance in the investigation of other three-flavor effects. 
To demonstrate this, we have plotted in Fig. [3] the sen- 
sitivity of the realistic ATLAS scenraio and several other 
experiments to sin 2 2^13. As examples for reactor ex- 
periments, we show Double CHOOZ \a, a nd a more ad- 
vanced setup with a 200 t far detector [46| . On the accel- 
erator side, we show simulation results for MINOS p7| . 
T2K [J, and NOi/A The sensitivity to sin 2 26 13 

is defined as the limit, which a specific experiment can 
set to sin 2 2013, assuming the true value is zero. The 
left edges of the bars correspond to the hypothetical case 
that only statistical errors are present in the experiment. 
The blue, green, and yellow bars reflect the limitations 
of the sensitivity due to systematical uncertainties, para- 
meter correlations, and the sign(Am 2 1 ) degeneracy 2 We 
can read off from Fig. [3] that atmospheric neutrino ex- 
periments can in general not compete with reactor and 
accelerator setups. Only Hyper-Kamiokande and the ILC 
scenario could improve the CHOOZ bound, though not as 



2 The octant degeneracy [45I is irrelevant here since we have as- 
sumed = 7r /^! an d the intrinsic degeneracy [49l | appears 
only in high statistics experiments such as a neutrino factory. 
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Figure 3: Sensitivity of various neutrino oscillation experi- 
ments to sin 2 2#i3 . The colored bars reflect the limitations of 
the sensitivity due to systematical errors, parameter correla- 
tions, and degeneracies. The right edges of the bars corre- 
spond to the sensitivity which is expected in reality. 



much as dedicated experiments. It is interesting to note 
that the ILC scenario has a significantly better sensitiv- 
ity than Super-Kamiokande, although it has only ~ 1, 350 
contained events and ~ 1, 720 upward going muons, 
while Super-Kamiokande has detected ~ 4,500 fully con- 
tained /i-likc events and ~ 2, 250 upward going muons, 
plus a large number of e-like events and partially con- 
tained events [HI]. This comparison confirms the state- 
ment that, to a certain degree, excellent resolutions can 
compensate for low statistics in atmospheric neutrino ex- 
periments [l2l ]. 



6. CONCLUSIONS 

We have shown that atmospheric neutrino interactions 
are an interesting by-product of operating the ATLAS 
detector. Although reconstruction of these events is re- 
stricted to phases where the LHC is not running, or only 
running at low luminosity, the number of events, in con- 
junction with the expected good energy and angular res- 
olutions, is sufficient to confirm atmospheric neutrino os- 
cillations and to measure the leading oscillation para- 
meters 023 and Ara^. However, the precision of these 




Muons 



Figure 4: Geometry of upward going muons. The white rect- 
angle schematically shows the ATLAS experiment as seen from 
a direction orthogonal to the detetor axis and to the muon 
direction. The shaded region represents the effective target 
volume V. 



measurements will probably not be competitive to that 
of existing and upcoming dedicated experiments. Under 
more optimistic assumptions, which might be realized at 
future ILC detectors, the errors on #23 and Anijj can 
become comparable to those of T2K, and, if #13 is very 
large, the sign(Am| 1 ) degeneracy might be resolved at 
the 2(j level. 

To study the reconstruction capabilities and sensitiv- 
ities of ATLAS in more detail, it will be necessary to 
perform detailed detector Monte Carlo simulations. We 
believe this to be definitely worthwhile, since a measure- 
ment of the atmospheric neutrino oscillation parameters 
with the technology of ATLAS would be a very inter- 
esting result in itself, even if the final error bars should 
be larger than those of other experiments. Moreover, the 
exploration of neutrino interactions in ATLAS would pro- 
vide valuable experience that could be of great interest 
in the design of future collider experiments and neutrino 
detectors. 
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Appendix A: GEOMETRY OF UPWARD GOING 
MUONS IN ATLAS 

As mentioned in Sec. [H the target volume for upward 
going muon events depends on the energy and on the 
zenith and azimuthal angles. It can be calculated from 
geometrical arguments as follows. We first integrate the 
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inverse of the Bethe-Bloch formula [5(| H3 to obtain the 
"muon range" L M , i.e. the average distance that a muon 
with energy can travel in rock. Wc then use Fig. [4] 
to calculate the effective target volume V in which a 
neutrino must interact to induce an upward going muon 
event in the detector. The white rectangle in the diagram 
represents the cylindrical ATLAS detector as seen from a 
direction orthogonal to the detector axis and to the muon 
direction, so that the angle T] between these two vectors 
lies in the drawing plane. The effective target volume, 
shown by the shaded region in Fig. dj is given by 

V = smr]\ ■ £ ■ 2r + L M | COS77I ■ nr 2 . (Al) 

We can calculate r\ from the azimuthal angle 4> and the 
zenith angle 9 by 

cos r\ = sin 9 cos (j>. (A2) 
Averaging over <j) yields 

V = ^- [ dcj)(2r£Jl- sin 2 9 cos 2 <j> + nr 2 sin 9 cos 4> \ ) 
27r Jo v 
ArPT 

= ££(sin0) + 2r% 1 |sin0|, (A3) 



with E(-) being the complete elliptic integral of the sec- 
ond kind [52| . 

When evaluating this expression, we assume a hori- 
zontally zylindrical geometry for ATLAS, with a length 
of 42 m and a radius of 11 m. For the geometrical ar- 
guments, the muon direction is taken to be identical to 
that of the primary neutrino, which is justified by the 
fact that upward going muons typically have very high 
energies 3> 1 GeV. The muon energy is calculated from 
the neutrino energy according to the empirical formula 



^ = K ^5.255-1.8191og 10 (fWMcV) + Q.3298^ , (A4) 



which we have obtained by fitting E^/E^ in a sample of 
1,000 upward going muon events in Super-Kamiokande, 
simulated with the NUANCE event generator [53j . 
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